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Abstract. — A family of random models for bosonic quasi-particle excitations, e.g. the 
vibrations of a disordered solid, is introduced. The generator of the linearized phase space 
dynamics of these models is the sum of a deterministic and a random part. The former 
may describe any model of identical phonon bands, while the latter is a c/-dimensional 
generalization of the random matrix model of Lueck, Sommers, and Zirnbauer (LSZ). 
The models are constructed so as to exclude the unphysical occurrence of runaway solu- 
tions. By using the Efetov-Wegner supersymmetry method in combination with the new 
technique of superbosonization, the disordered boson model is cast in the form of a su- 
permatrix field theory. A self-consistent approximation of mean-field type arises from 
treating the field theory as a variational problem. The resulting scheme, referred to as a 
coherent potential approximation, becomes exact for large values of A^. In the random- 
matrix limit, agreement with the results of LSZ is found. The self-consistency equation 
for the full ^/-dimensional problem is solved numerically. 



1. Introduction 

Small oscillations about the stable equilibrium of a many-body ground state are 
quantized as bosonic quasi-particles or bosons. In various physical contexts the lin- 
earized equations of motion for such excitations are known as the random phase ap- 
proximation, or RPA equations for short [21. Concrete examples are furnished by the 
vibrational excitations of a solid, the spin waves of a magnet, the electromagnetic 
modes of an optical medium, or the density oscillations of a Bose-Einstein condensate. 

Constrained by the requirement of dynamical stability, the Hamiltonian H of any 
vibrational or quasi-boson system of the mentioned kind must lie in a positive cone, 
S", of so-called elliptic symplectic generators. It should be stressed that although H 
is Hermitian as an operator in Fock space, the quantum- to-classical mapping sends H 
to an RPA generator X which is in general neither Hermitian nor anti-Hermitian as 
a linear operator on the classical phase space. In view of this, a distinctive feature 
of the set of elliptic generators X E tS' h that they can be brought to diagonal form 
(with real frequencies, corresponding to stable oscillatory motion) by real Bogoliubov 
transformations, i.e., by conjugating with elements of the real symplectic group Sp^ . 
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In this general setting, our goal is to investigate what happens with observables 
such as the spectral statistics and the transport properties when the bosonic system is 
strongly disordered. In particular, we wish to understand whether there exists some 
low-energy universality, possibly of an unusual type, due to the interplay between dis- 
order and the geometry of the cone ^ . (For example, a high degree of low-temperature 
universality is known to be observed @ in strongly disordered solids as well as amor- 
phous or glassy systems.) Motivated by this question, in the present paper we initiate 
the study of a class of semi-realistic random matrix models for disordered bosons. 

By construction, the probability measures of the models we propose are supported 
on S . Thus, unlike [3], the unphysical behavior of runaway motion associated with 
complex frequencies is excluded. Our RPA generators X = K^R have the particular 
feature of being sums of a deterministic and a random part. For simplicity we focus 
here on models without time-reversal symmetry, although TR-invariant models of a 
similar kind can be treated with little extra effort. By using a variant of the Efetov- 
Wegner supersymmetry method, we will derive an approximation for the density of 
states of mean-field or self-consistent type, reminiscent of the so-called 'coherent po- 
tential approximation' (CR^) [|4l|5l|6l. 

In the longer term, the goal is to develop a description of our disordered boson 
models by field theories of the non-linear sigma model type. (As is well known [7J, 
such a description has proven very useful for the case of disordered fermions). In that 
formulation, universality (if any) is expected to emerge whenever the renormalization 
group flow gets attracted to a few -parameter manifold of renormalizable field theories. 

The plan of the paper is this. In Section [2] we outline the basic setting and in partic- 
ular, we review the notion of positive cone of elliptic symplectic generators. We also 
introduce the random models to be considered and give a summary of the analytical 
results obtained. Section |3] is concerned with the derivation of the coherent potential 
approximation for our models. For pedagogical reasons, we first discuss the zero- 
dimensional case in some detail. The extension to d dimensions is given in Section IH 
There we also show some results for the numerical solution of the CR^ equation. 

2. Setting, model, and results 

In this section, we start with some background on mathematical formulation and de- 
scribe a class of random models which are tractable by the superbosonization variant of 
the Efetov-Wegner supersymmetry method. We then give a summary of our analytical 
results, relegating the presentation of numerical results to the end of the paper. 

2.1. Setting. — Let {aj,a^}yter be a set of creation and annihilation operators for 
bosonic quasi-particles with quantum numbers ^ G P. (For example, F might be a 
discrete set of momenta selected by periodic boundary conditions in a finite box.) 
Such operators span a Hermitian symplectic vector space W := span£{a^,aj}^£r with 
symplectic form [ , ] : x V7 — )■ C defined by the canonical commutation relations 
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We now assume that we are given a linear Hamiltonian dynamics on W. This may be 
interpreted either as a linear Hamiltonian flow on W viewed as a classical phase space, 
or as a quantum time evolution on W viewed as a subspace of the associative algebra 
of polynomials in a^, a J (the so-called Weyl algebra). In either case, the dynamical 
equations are 

^ = L i'^l'Yk'k + av^k'k) , ^ = L {'^I'^k'k + cik'yk'k) ■ (2.1) 

(The bar means complex conjugation.) In order for the canonical commutation rela- 
tions to be invariant under the dynamics, we require that 7^^/ = — F/t'/t ^^id Z^^/ = Z^/^. 
Thus Yj^j^i are the matrix elements of an anti-Hermitian matrix Y = —Y\ while Z^^yf are 
those of a complex symmetric matrix Z = Z^. Altogether, these conditions mean that 

X := f) (2.2) 

is the generator of a symplectic transformation. More precisely, defining the Lie alge- 
bra, sp, of the complex symplectic group by the linear condition 

X = -JX'J\ -^=(-1 o)' ^^-^^ 

X lies in a non-compact real form sp^ C sp determined by 

X = -E3X%, S3=(J (2.4) 

It should be mentioned that this description is appropriate in the absence of time- 
reversal invariance. If time reversal is a symmetry of the physical system, then the 
time-evolution generator X is subject to additional complex anti-linear conditions. 

Equations (|2.1I) arise as the equations of motion for a system of non-interacting 
bosons with Hamiltonian 

H = iJ'^Y,k,k' i^kk^alcik' + {^kk'Cikak' - \ZkMalal) (2.5) 

and dynamics ihd= [H,a\. Alternatively, one may imagine that they arise as an ap- 
proximation to the collective motion of an interacting many-particle system; as a par- 
ticular example we mention density oscillations of a fluid. In the latter case, equations 
(12.11) go under the name of random phase approximation (RPA). 

The characteristic frequencies of the dynamical system (12.11) or equivalently, the 
single-boson energies of the Hamiltonian H, can be computed as the eigenvalues of 
X. Owing to the symplectic condition X = —JX^J^ the characteristic polynomial 
satisfies Det(A —X)= Det(A +X). The eigenvalues of X therefore come as pairs ±A. 

If X lies at some random position in the real symplectic Lie algebra 5pR , then its 
eigenvalues will typically be complex, since X E sp^ is neither Hermitian nor anti- 
Hermitian. In the present context, however, complex eigenvalues are forbidden, as 
they correspond to the unphysical situation of unstable motion. In fact, the physical 
requirement of stability of the RPA dynamics dictates that the spectrum of X must lie 
on the imaginary axis, so that the normal modes of the bosonic system are vectors in 
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W with periodic time dependence (o« e^^^'). Moreover, the second-quantized Hamil- 
tonian H in (12.51) must have a ground state in Fock space. By standard considerations 
of linear algebra, all these stability conditions are met if and only if X lies in the set 

^■={Xesp\ iZ^X > 0}. (2.6) 

We refer to S' as the positive cone of elUptic generators in sp^ . It is a fact that every 
X G (f can be brought to diagonal form by a real Bogoliubov transformation, i.e. an 
element g of the real symplectic group Sp^ , which is defined by the condition 

J{g-'p-'=g = Mg-')^^3. 



2.2. The model. — In the present paper we consider RPA generators X, or equiva- 
lently Hamiltonians H, which are a sum of two parts: 

X = K + R. (2.7) 

The term K is the deterministic (i.e., non-random) part of X. While the formalism 
developed below can in principle handle any choice of K, the explicit calculations 
presented in Section l4n will be carried out for a simple concrete model of K with unit 
mass matrix and elastic constants given by a discrete Laplacian. A precise description 
of the concrete model for K is as follows. 

2.2.1. Deterministic part. — Let A = Z"^ be a cubic lattice in d space dimensions 
and associate with each site j E A the operators and aj for boson creation and 
annihilation, respectively. We then take the second-quantized Hamiltonian to be 

H = hvY,ci]aj-^ £ {aj+a]){aj,+a],) (2.8) 
jeA UJ') 

where the sum for the second term on the right-hand side is over nearest neighbor pairs 
of sites J, / of A. Such a Hamiltonian is easily diagonalized by Fourier transforming 
to momentum space. The spectrum of single-boson energies e{k) as a function of the 
wave vector k= {k\_,...,kti) comes out to be 

1 

eik)=Tiv^/l-Ak, Ak = -Y,cos{k). (2.9) 

" 1=1 



Note that £{k) ~ hv\k\/y/2d for small \k\ = ^kj + ... + kj, which tells us that the 

speed of sound in units of the lattice spacing is v/ \/2d. 

By computing the RPA generator from the dynamical equation iha = [H, a] for 
a = a^j and a = aj we obtain the expression 

^f':^ .-^). (2.10) 



2 \ A -2 + A 

where A is the scaled lattice Laplacian (with diagonal part removed) which has eigen- 
value spectrum A^^. Next, we tensor up the model by introducing A'^ identical bands. 
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Mathematically speaking, we pass from the symplectic vector space (for each j G A) 

spanc{a^-,aJ}~C2 

to the tensor product Wj := g) C^^ and take the generator K to be 

K■=K^^ldcN. (2.11) 

This means that creation operators a j ^ and annihilation operators „ get an extra band 
index n = 1, . . . , A^. Note that in the physical setting of lattice vibrations a reasonable 
choice of in J dimensions is N = d due to the vector nature of lattice displacements. 

2.2.2. Random part. — We turn to R, the second term in (12.71) . which is random. A 
particular feature of our disordered model is that we take R to be diagonal in the sites 
J G A of the lattice. For simplicity we begin the discussion with the very special case 
of A consisting of just a single site. The full model to be discussed later is obtained by 
repeating the single-site discussion at every site of A = Z"^. 

With the single site of the lattice we associate a Hermitian vector space W = C^^ 

with symplectic structure 7 = ^ ^ . In order to implement the positivity con- 

dition [see Eq. (12.61) 1 for X to be in the cone <f , we let 

R = -iT.^L^L, S3=(o^ (2.12) 

where L is a rectangular linear operator 

L:W^V, y = C^, (2.13) 

mapping W into an auxiliary vector space V . The dimension M is a parameter of our 
model. It may be bigger or smaller than IN. A special role is played by the choice 
M = IN, as this is the minimal dimension for the operator R to have full rank. 

It is easy to see that for R = —iLjL^L the symplectic condition R — —JR^J^ holds 
if and only if L satisfies the reality condition 

^ = "'' o)- <2.14) 

This condition fixes a real form, say U^, of the complex vector space U = Hom(W^, V). 
Note that if L'f'L has full rank then iZ^R = L^L>0 and RES'. 

Disorder is introduced by declaring the matrix elements of L to be Gaussian random 
variables. More precisely, we define the probability measure for L G t/R as 

d^{L) =Ce-^^'^^^dL, (2.15) 

where dL is Lebesgue measure on the normed vector space [/r and C is a normalization 
constant. The parameter bis a measure of the disorder strength. We mention in passing 
that the model forX =R = -iZ^L^L with probability measure (12.151) (and M > 2N) is 
equivalent to the random matrix model studied in BHl by different methods. 

Finally, we describe the generalization to an arbitrary lattice or graph A. With each 
lattice site j G A we associate one copy Wj of the Hermitian symplectic vector space 
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C^^. The total physical space then is the orthogonal sum W = (BjeAWj . Note that 
dim(V7) = 2A'^|A| where |A| denotes the number of sites of A. The full generator of 
the dynamics is X = K + R where the deterministic part K may in principle be any 
element of the positive cone £' (W) . For concreteness we let A = and take K to be 
the generator described in Section I2.2.1I The random part i? is a sum R = J^j Rj of 

Rj = —iL^LjLj made from independent and identically distributed random operators 
Lj . In other words, the distribution for R is given by the product distribution 

d^A{L) = lld^{Lj). (2.16) 

2.3. Statement of result. — While our interest will ultimately be in correlation func- 
tions and transport properties, we here take a first step by studying the average resol- 
vent of the time-evolution generator X : 

g{z) = (2^|A|)-iE(Tr(z-X)-i) , (2.17) 

where the symbol E(. . .) means the expectation value with respect to the probability 
measure (12.161) . Notice that by the symplectic property X = —JX^J^ the resolvent 
satisfies Tr(z — X)^^ = Tr(z4-^)^^ so g{z) = —g{—z) is an odd function of the fre- 
quency parameter z ■ Because the support of our probability measure is contained in the 
positive cone of elliptic elements, <f, the eigenvalue spectrum of the random operator 
X is always imaginary and g{z) is analytic in the right and left halves of the complex 
z-plane. In the following we assume £Hez > 0. It is a standard fact that the local density 
function p of the characteristic boson frequencies ft) can be computed from 

p{a>) = n^^ lim 9^e^(±ioj + e). 

We now come to our main result. Fixing the ratio 

a-.^MjW, (2.18) 

we take the large-A/^ limit of the model with dynamical generator X=^ + i?onA — Z"^ 
as described above. We then claim that in this limit g[z) is expressed by 

[0,2k]'' 

where the complex and energy-dependent quantity p is a solution of the self-consistency 
equation 

1-" f ''"'^ P + v(l-iA,) ^^^^^ 



b p J {InY z^+p^+pv{2-Ak)+y2{l-Ak)' 

[0,2;r]'' 

p plays the role of a 'self energy' or 'coherent potential'. 

We briefly discuss some features of the solution in two extreme cases. There is only 
one relevant parameter, b/v. In the limit of weak disorder (b — )■ 0) one infers that 
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p — )■ and 

S{Z) = T^ / . ..1. . . (2.21) 



{2kY J z^ + v\\-^k) 

is simply the Cauchy transform of the local density of frequencies of the deterministic 
generator K. On the other hand, for strong disorder {b — > oo) the coherent potential 
p becomes large and the system (|2.19[ 12.201) simplifies to 

z^ + p^ b p z^ + p^ 
A special situation arises for a = 1. In this case it follows by a short computation from 
(12.221) that the scaled function g{x) := ibg{ibx) satisfies an equation, 

-1 

which was derived and solved by Lueck, Sommers, and one of the authors |[8l. The 
analysis of |l8l shows that the density of states behaves as x^^^^ near x = in this case. 
Numerical results for I and for the more general situation of intermediate disorder 
strengths will be presented in Sections [331 and |4~T1 



3. Solving the model in zero dimension 

In the sequel we explain how to arrive at our main equations (12.191) and (12.201) . For 
pedagogical reasons we describe the method first for the simple situation of a single 
site (the zero-dimensional case). Throughout this section we let W = and V = C^. 

3.1. Resolvent as Gaussian integral. — Our plan is to compute the average trace of 
resolvent (12.171) by a variant of the Efetov-Wegner supersymmetry method. The first 
step of this method is to express Tr(z — X)^^ as a Gaussian Berezin (super-)integral. 
To get started, we use the elementary identity 



Tr(z-X)-i = — 



Det(zi-X) 

— (3 I) 

Det(zo-X)' ^"-'^ 



Zl=ZO = 

and then write each of the two determinants as a Gaussian integral - using ordinary 
integration variables for the determinant in the denominator and anti-commuting vari- 
ables for that in the numerator. 

In the case of the ordinary Gaussian integral, there exists a convergence issue be- 
cause the elements X G spR have indefinite real and imaginary parts in general. It is 
therefore crucial that all our generators X, constrained to lie in the positive cone S', sat- 
isfy the inequality iL^X > 0. To take advantage of this positivity property, we express 
the determinant as follows: 

Det(z-X) =Det(z-i^ + iE3LtL) =z2^-^A(z), (3.2) 
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Note that owing to £Hez > and iL^K > 0, the 2 x 2 matrix of operators in (|3.3I) has 
positive real part. 

Next, we introduce symmetric complex bilinear inner products ( , ) for each of the 
two vector spaces V and W. These inner products are consistent with the Hermitian 
structures of V and W in the sense that, e.g. for V, the sesqui-bilinear form (v, v') t-)- 
(v,v') agrees with the Hermitian scalar product of V. We then express the reciprocal 
determinant 1/A(z) as an integral over two complex vectors v eV and w gW : 

A(z)^^ = y"e^^(^'"^)+'^"'(^3vi',H'z-A:w)-i(v,Lw)-i(H',L+v) q ^-j 

where it is understood that we are integrating with the product of Lebesgue measures 
for V and W. The normalization is chosen in such a way that / e^(^''^')^(''^''"') = 1. We 
emphasize that the integral (13.41) converges absolutely due to Rez > and K E S'. 

In the case of the determinant itself we integrate in the sense of Berezin (i.e., we 
actually differentiate) with respect to two independent vectors /3 and 7 whose compo- 
nents are Grassmann variables: 

= yez(j6,/3)-iz-i(Z3r,7z-^^r)+i()3,ir)+i(y,i"^/3)^ (3_5) 

Again, it is understood that we are integrating with the flat Berezin form, i.e., the 
product of all partial derivatives w.r.t. the Grassmann variables. The bar in the present 
instance means nothing but independence, e.g., of j8 from /3. 

We now multiply the two Gaussian integral formulas (13.41) and (13.51) and take the 
disorder average inside the absolutely convergent integral to obtain 

\Mzo)J J 

X E j^e"i(v,M-i('^',-'^"' v)+i(j8,Lr)+i(y,Llj8) j _ ^2.6) 

This integral representation is a suitable starting point for further analysis. 

3.2. Taking the disorder average. — Next, we compute the disorder expectation 
value in (|3.6I) . For that we introduce the quadratic quantities 

2 := v(w, •) +v(Eiw, •) +/3(r, ■) -i8(Eir, ■), 
Q* := w (v, ■) +EiM) (v, ■) + r(i8, ■) - Eif(/3, ■), 

where Q is meant as a linear transformation from to V with coefficients in the even 
part of a Grassmann algebra, and similar for Q* with the roles of W and V reversed. 
We then have 

^ j^g-i(v,Lw)-i(w,L' v)+i(i3,Lr)+i(r^L"'"j8) j ^ J q-^TH^'^' Q+Q*L) ^^(^^ ^ Q-^TrQ*Q^ 

by completing the square and shifting variables. 

At this point we make the observation that Tr Q* Q depends on v, v, j8 , j8 only through 
scalar products such as (v,v), (v, v), (v,/3). These share the feature of invariance under 
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(3.7) 



the group Om of real orthogonal transformations of V = C^. It will be useful to 
organize all these OM-scalars into a supermatrix: 

/(v-,v) (v-v) (v-i8) -(v-J)\ 

(v,v) (v,v) (v,i8) -(vj) ^fPoo Poi 

(j8,v) (j8,v-) (i8,/3) -\Pw Pn 

\{P,v) (j8,v-) -{I5J)J 

Two of the matrix entries vanish since (j8 , j8 ) = — (j8 , /3 ) = and, similarly, (j8 , j8 ) = 0. 
We also have (j8, j8) = — (/3, j8). We further note the expression 

STr P = TrPoo - TrPi i = 2(v, v) - 2(j8 , j8 ) 

for the supertrace of P. 

There exist certain linear dependencies amongst the matrix elements of P. To de- 
scribe them we need the operation Ph^ P^^ of taking the supertranspose: 

/lO ^11 

With its help we can formulate the symmetries of P as follows 



P^ 

^00 
M)l 



^ 10 
Ml 



(7 = diag(c7i,ic72), (7i 



-i 

1 



We now arrange the remaining integration variables w, w, y, and 7 in the form of 
rectangular supermatrices: 

V(^ir,-)/ 

More precisely, ^ is to be viewed as a linear mapping from W into the superspace C^l^ 
with Grassmann-even resp. Grassmann-odd matrix coefficients on the even resp. odd 
positions of Hom(W, C^'^). It is easy to check the identity 

Tr2*2 = TrW^, 

which lets us re-express our disorder average as 



E e 



i(v,Lw) 



,L■'v)+i(i3,Lr)+i(y,Lt]3)^ ^ -^TrW^P 



(3.8) 



3.3. Eliminating ^, ^. — The next step is to carry out the integral over w, w, 7, and 
7, thereby eliminating ^ and from the calculation. This will be straightforward to 

do because the dependence on these variables is Gaussian. 

As a preparatory step, we verify from K = —JK^J^ the relation 

i{'L3W,w-Zo'Kw) - i(L3 7, 7- zr'^r) = ^TriiZ^^z^'i' -i^sK^r'^V), (3.9) 
where 

£ = diag(zo,zo,zi,zi) , '^3 = diag(l, -1, 1, -1) . (3.10) 
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By using the results (|3.8[|3.9I) in equation (13.61) we then arrive at our next formula: 

e( ^^^^ = f Q-jSTT{zP)+^Tr{iI.3^T3'i'-iL3K^z-^'i'-^^P'i') 

The integral on the right-hand side is still over the original variables v, v, /3, j8 in P 
and w, w, y, f in ^. Finally, by using a standard formula for Gaussian Berezin 
superintegrals we perform the integral over ^ and ^. This results in 



E 



The superdeterminant here is over the tensor product space W (g) C^'^. We recall that 
the superdeterminant of a supermatrix is defined by 

/A B\ Det(A) Det(A-5D iC) 

SDet ' ' - V ^ _ \ ^ 



C DJ Det(D-CA I5) Det(D) 

The integral above is still over the variables v, v, /3, and j8 entering via their scalar 
products into the supermatrix P. By scaling these integration variables so that P — )• 
2NP I b, we obtain the following expression for the generating function of our problem: 

X J e-^^^''^''^^SDet-^/^{-iL3®zT3+i^3K®ld2\2+ldw^zP). (3.11) 
The symbol Id2|2 stands for the identity in superspace C^'^. 

3.4. Reduction by superbosonization. — Superbosonization is a change of vari- 
ables \Qj which lets us switch from integrating over a large number of vector-type 
variables, to integrating over a smaller number of matrix-type variables. In the present 
context these are the components of the vector variables v, v, /3 , j8 and the matrix ele- 
ments of the supermatrix P, respectively. Such a reduction of the number of integration 
variables does not come for free but requires the integrand to be invariant under one of 
the Lie groups GL, O, or Sp. There exists a version of superbosonization for each of 
these classical Lie symmetries. As we have seen, our integrand is expressed in terms 
of quadratic invariants of the orthogonal group Om • Therefore we now recall from [|9l 
the superbosonization identity for the case of O^-symmetry. 
The OM-superbosonization identity reads 

J F{P{v,v,l5,~P)) = J D^{P) SDet^/2(P)F(P) , (3.12) 

where on the left-hand side we integrate with the flat Berezin form 

I I dVmdVm — 
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and on the right-hand side the Berezin integration form is 

D^{P) =DPSDet^/^{P), (3.13) 

where DP is still the flat Berezin form (i.e., the product of differentials for the even 
variables and partial derivatives for the odd variables). The domain of integration for 
the so-called boson-boson block Pqq [see Eq. (13.71) 1 is the space of positive Hermitian 
2x2 matrices Pqq = Q subject to Q — OiQ^Oi . In the fermion-fermion sector, the inte- 
gration domain is the space of unitary 2x2 matrices Pn =U subject to the symmetry 
relation U = 02U^02 ■ These matrix spaces are diffeomorphic to the symmetric spaces 
GL2(M) /O2 and U2/USP2 respectively. 

By applying the superbosonization identity (13.121) to the integral representation (13.111) . 
we obtain our final result for the generating function: 

a{z) = j D^{P) e-^STrP sDet^/2(/') 

X SDet"^/2(-iE3®zT3 + iE3^®Id2|2 + Mvv®P). (3.14) 

Notice that a substitution P ^ z^P was also made. By the relation SDet^/^(|^^/') = 
izi/zo)^SDet^/^{P) this removes the multiplicative constant (zo/zi)^ from (13.111) . 

The result (13.141) is exact and mathematically rigorous for M >2. (In the present 
case of OM-symmetry the superbosonization identity fails for M = 1; see [i9|.) From it 
we get the average trace of resolvent by differentiating at coinciding points zo = zi'- 



E(Tr(z-X)-i) = |- 



a(diag(zo,zo,2i,2i))- 



3.5. Random-matrix limit. — To conclude this section we consider the special limit 
of vanishing deterministic generator ^ = 0. In that case our expression simplifies to 

a{z) = J D^{P) e-^STrP sDet^/\P) 

X SDet-^/2 _ ji^^) SDet-^/2 (P + izTs), (3.15) 

where all superdeterminants and supertraces are over C^'^. Recalling the parameter 
a = M/2N we see that our integral is of the form 



aiz) = jo^iP) e- 



-NF(P) 



with 



F (P) = Z?" ^ STr P - a In SDet (P) + i In SDet (P - izT3 ) + 5 In SDet (P + izTi) . 

We now investigate the random-matrix limit — )• 0° with a = M/2N held fixed. In this 
limit the integral for Q.{z) can be computed by the saddle-point or Laplace method. 
By the principles of supersymmetry, the leading contributions to the integral at zo — Z\ 
can be shown filOJ to come from saddle points which are multiples P = pld2|2 of the 
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2 4 1 



U> UJ 

Figure 1 . Density of eigenfrequencies in the random-matrix limit K = 0. 
The parameter values are ft = 1 and a = 2 (left), a = 0.75 (right). In the latter 
case there exists a Dirac-delta peak (1 — a)5{(o) due to zero modes. 



identity. We here omit the details of the calculation and present only the outcome. By 
execution of the saddle-point method we find that 

g{z) := lim (2A^)-iE (Tr(z-X)-i) = (3.16) 

where p is a solution of the saddle-point equation 

1 a 1/2 1/2 



b p p~iz p + iz 
In Section [23] [see Eq. (12.221) 1 we already remarked that for a = 1 this is equivalent to 
an equation analyzed and solved in [8J. Hence in what follows we focus on a 7^ 1. 

We begin with the case a > 1 . A plot of the density of states for a = 2 is shown in 
Figure [U We see that there is a gap at low frequencies. This feature can be understood 
in the same way as the Marcenko-Pastur law [ill II for rectangular Wishart matrices. 
Indeed, recall that our random generator is R = —iL^L^L where L E Hom(W^, V) is 
rectangular of size M x 2N. The non-zero eigenvalues of R coincide with those of the 
operator —iLL^L^ but the latter has M — IN = IN {a — 1 ) additional eigenvalues at zero 
by rank-nullity. In the large-A^ limit the level repulsion due to this macroscopic number 
of zero modes produces a spectral gap of size proportional to {a — I). 

The gap closes as a approaches unity, leading at a = 1 to the situation investigated in 
BHl . For a < 1 it is the operator L^L which by rank-nullity has 2N — M — 2N{ I— a) zero 
modes, and the same goes fori? = —iLt,L^L. Therefore the density of states contains a 
Dirac mass (1 — at zero in this case. A plot of the density of states for a < 1 is 

shown in Figured! where see that the DOS approaches a finite value at zero frequency. 
(The argument of macroscopic level repulsion does not apply here, as the operator R is 
neither Hermitian nor anti-Hermitian.) 

The discussion above is concerned with the so-called bulk scaling limit. Another 
limit of interest is the edge-scaling limit at z = where one sends — t- 0° while keeping 
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z\/N fixed. For a = \ this limit was thoroughly investigated in BHl, while for a > 1 the 
situation is trivial because of the absence of states at z = 0. 

For a < I the edge-scaling limit was studied in iflOl . The Hessian of the function 
F{P) at the saddle point P = pld2\2 has eigenvalues of order z. Therefore, in edge 
scaling z ~ N^^^^ this saddle point is not isolated and one has to work with a whole 
supermanifold of saddle points. (Technically speaking, the saddle-point supermanifold 
is a Riemannian symmetric superspace OSp/GL of type CI|DIII.) The law for the 
density of states in the limit N ^ oo turns out [jlOl to be the universal law for systems 
of class D in the symmetry classification of [Il2l . 



4. Going beyond zero dimension 

We now turn to the <i-dimensional model described in Section 12.21 The procedure 
of deriving the coherent potential approximation (I2.19[ 12.201) for this model remains 
essentially the same as before. Again, our first step is to express the determinants in 
(13.11) as Gaussian integrals over vector variables v, v, /3, j8 for V = (Bj^A^j and w, w, y, 
fforW = (BjeA^j ■ The Gaussian integral representation has the effect of factorizing 
the independent random variables associated with different sites of the lattice. The 
disorder average can therefore be carried out for each site separately. By the local 
Om gauge symmetry of the model, the integrand after disorder averaging depends only 
on Om gauge invariant combinations of the fundamental variables v, v, /3, j8. These 
organize into supermatrices P as before. Thus we introduce such a supermatrix Pj for 
each site j E A and switch to integrating over Pj by the superbosonization formula 
(13.121) . Because the dependence on the variables w, w, 7, 7 is still Gaussian, they can 
again be integrated out to produce a superdeterminant. In this way we obtain 

a(i) = [ ]\Dii{Pj) e-^STrP, sDet^/2(/',) 

X SDet-^/2 (^io^Kn^i ®Id2|2 + £/n^- ^Z'; -icTs ^Its)) , (4.1) 

where denotes the orthogonal projector from W/v=i onto (Wj)^^!- The only dif- 
ference of any essence from our earlier result (13.141) is that the integral now is over 
a field of supermatrices {Pj} j^a instead of a single supermatrix P. The operator 

Y^ji^j^Pj — i<73 ^zTi,) is diagonal on W but (for generic Pj) non-diagonal on C^'^. 
On the other hand, the operator iOsKi ® Id2|2 is trivial on superspace but couples the 
sites of the graph A. The inverse square root of SDet is raised to the A^* power because 
each of the bands of the deterministic limit contribute the same factor. 

We now face the task of analyzing the model (|4.1I) by the field-theoretic methods of 
gradient expansion and renormalization. (Note that a closely related problem has al- 
ready been tackled in [HOl .) Hoping to make progress with this in a future publication, 
we here take a first step by computing the local density of states. 
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Figure 2. DOS for J = 1, a = 0.75, v = 1 and = 0.63 (left), b = 0A5 (right). 

4.1. DOS for the concrete model. — Let us finally work out the mean-field solution 
of the model (14.11) with deterministic generator as defined in (12.101) . Writing the 
integrand as e^^^ (for M = IN a) we take the general variation of F: 

dF^b-^ Y^. STr 5Pj - a STr i dPj 

+ \ £ . STr (n,- ® dPj) {iOsKi ® 1 + £^ ( 1 ® - i(73 ® f T3) ) " ^ 

For large N we expect the field integral to be essentially given by a spatially homoge- 
neous saddle point Pj = pld2\2 (independent of j E A) and small fluctuations around 
it. Therefore, after setting zo = zi = z v^/e look for solutions of dF = of this very 
form. The variational equation dF = then reduces to an equation of self-consistent 
mean-field type: 

0=1-^ + 1 /■-^LTrl^"'^ + ^ + ^"i^^^ "^V^^ 

b p i] {InY V -\^^k iz + p + V-\v^k) ' 

where we have used the property that the Laplacian A is diagonal with eigenvalues A^ 
in momentum space. By evaluating the trace of the matrix inverse we immediately 
arrive at equation (12.201) . Within this mean-field (or coherent potential) approximation 
scheme, we obtain the expression (12.191 ) for the average resolvent trace g{z)- 

Let us finish by showing some numerical results for the model in dimension d — I. 
In this case the density of states from (|2.21l) for the pure system {b — 0) is 

As is seen in Figure[2l the van Hove singularity at co = y/lv is still visible for ^ = 0. 15 
(and a = 0.75). As the disorder strength b is increased, the bulk of the spectrum is 
pushed to higher frequencies and a peak begins to develop at small frequencies (see 
the plot for b = 0.63). At values of b much larger than the sound velocity v we recover 
the random-matrix limit shown in Figure [TJ 
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